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Abstract 

We analyze perturbed version of the complex Toda chain (CTC) in an at- 
tempt to describe the adiabatic A'^-soliton train interactions of the perturbed 
nonlinear Schrodinger equation (NLS). We study perturbations with weak quadratic 
and periodic external potentials by both analytical and numerical means. The 
perturbed CTC adequately models the A'^-soliton train dynamics for both types 
of potentials. As an application of the developed theory we consider the dynam- 
ics of a train of matter - wave solitons confined in a parabolic trap and optical 
lattice, as well as tilted periodic potentials. 

1 Introduction 

The A^-soliton train interactions for the nonlinear Schrodinger equation (NLS) and 
its perturbed versions 

iut + ]^Uxx + \u\'^u{x,t) =i€R[u], (1) 

started with the pioneering paper T , by now has been extensively studied (see |21 El 
01 13 12] and references therein) . Several other nonlinear evolution equations (NLEE) 
were also studied, among them the modified NLS equation El El El i some 
higher NLS equations the Ablowitz-Ladik system T* and others. 

Below we concentrate on the perturbed NLS eq. 1^. By A^-soliton train we mean 
a solution of the (perturbed) NLS fixed up by the initial condition: 

u{x,t = 0) = Y,ul^{x,t = Q), 4«(x,i) = ^^^, (2) 

^ cosh 2fc 
k—l 

Zk{x,t) ^2vk{x~ik{t)), ik{t)^2^ikt + £,k,o (3) 

4>k{x,t) = ^Zfe + Skit), 6kit) - Wkt + 4.0, (4) 

Each soliton has four parameters: amplitude Vk, velocity fik, center of mass position 
^fc and phase Sk- The adiabatic approximation uses as a small parameter eg ^ 1 the 
soliton overlap which falls off exponentially with the distance between the solitons. 
Then the soliton parameters must satisfy [J: 

\vk - Vol <t: vo, - Mo| < A*o, l^^fc - J^ol 16+1,0 - 6,o| > 1, (5) 

where i^o ^ jf 12k=i '^k, and — jj J2k=i l^k are the average amplitude and velocity 
respectively. In fact we have two different scales: 

I I 1/2 I I 1/2 1^ ^ \ -1 

- i^ol — £0 1 lMfe~Mo| — £0 ' |?/c+i,o - 4fe,o| — £0 • 



1 



One can expect that the approximation holds only for such times t for which the set 
of AN parameters of the soliton train satisfy ©. 

Equation Q finds a number of applications to nonlinear optics and for R[u\ =0 is 
integrable via the inverse scattering transform method |13l I14| . The 7V-soliton train 
dynamics in the adiabatic approximation is modelled by a complex generalization of 
the Toda chain [TS] : 

^ = 161.', {eQ^^^-Q^ - eQ^-Q^-^) , J = 1,...,N. (6) 
The complex-valued Qk are expressed through the soliton parameters by: 

Qfe(t) = 2iAo6(i) + 2fcln(2i/o)+*(fc7r-4(t)-<^o), (7) 

where i5o = and Ao = /io + *f^o- Besides we assume free-ends conditions, 

i.e., e~'^° = e'3«+i ee 0. 

Note that the A^-soliton train is not an A^-soliton solution. The spectral data of 
the corresponding Lax operator L is nontrivial also on the continuous spectrum of L. 
Therefore the analytical results from the soliton theory can not be applied. Besides 
we want to treat solitons moving with equal velocities and also to account for the 
effects of possible nonintegrable perturbations R[u]. 

The present paper extends the results of several previous ones: see Refs. 21 
01 El E] • Recently with the realization of Bose-Einstein condensation of dilute 
atomic gases it became important to study NLS equation with additional potential 
term iR[u] = V{x)u{x,t), see ^lEl- We continue the analysis in JBJ and in our 
more recent reports ^^1; of the corresponding perturbed CTC (PCTC) model for 
quadratic and periodic potentials V{x). Our results give additional confirmation of 
the stabilization properties of the periodic potentials observed in |20[I21| in a different 
physical setup. 

We also pay attention to the so-called tilted periodic potentials, which are super- 
positions of periodic and linear potentials. The effect of the linear potential is that it 
can, if it is strong enough, overcome the confining effect of the periodic potential. As 
a result we can have one or more of the solitons extracted out of the train. We also 
demonstrate that the PCTC provides an adequate description also for these types of 
potentials. 



2 The importance of the CTC model 

The fact [El E21 that the CTC, like the (real) Toda chain (RTC), is a completely 
integrable Hamiltonian system allows one to analyze analytically the asymptotic be- 
havior of the A^-soliton trains. However unlike the RTC, the CTC has richer variety 
of dynamical regimes 01211221 such as: 

• asymptotically free motion if Vj ^ Vk for j ^ k] this is the only dynamical 
regime possible for RTC; 

• iV-s bound state if ui = • • • = vn but (^k 7^ Cj for k ^ j; 

• various intermediate (mixed) regimes; e.g., if ui = W2 > ■ • • > I'A' but (^k ^ Cj 
for k j then we will have a bound state of the first two solitons while all the 
others will be asymptotically free; 

• singular and degenerate regimes if two or more of the eigenvalues of L become 
equal, e.g., Ci = C2 • • ■ and Q ^ Cfc for 2 < j ^ k. 

By Ck — + iwk above we have denoted the eigenvalues of the Lax matrix L in 
the Lax representation L-r ~ [Af, L] of the CTC where: 

N N-1 

L ^^^kk + ak{EkM+i + Ek+uk), (8) 

k=l k=l 
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and the matrices Ekp are defined by {Ekp)ij = SkiSpj. The eigenvalues of L are 
time independent and complex-valued along with the first components rjk = z^'''^ of 
the normalized eigenvectors of L: 

LZ^"^ = Ck^'^K (Z<'=), Z^")) = Skm. (9) 

The set of {Ck = Vk + iwk, Tjk = cTfc + iOk} may be viewed as the set of action-angle 
variables of the CTC. 

Using the CTC model one can determine the asymptotic regime of the A^-soliton 
train. Given the initial parameters /i/c(0), i'fe(O), ^^(O), (S/£(0) of the iV-soliton train 
one can calculate the matrix elements bk and afc of L at t = 0. Then solving the 
characteristic equation on L|t=o one can calculate the eigenvalues Cfe to determine 
the asymptotic regime of the A^-soliton train [21 E] . Another option is to impose on 
C,k a specific constraint, e.g. that all C,k be purely imaginary, i.e. all Vk = 0. This 
will provide a set of algebraic conditions L|t=Oj and on the initial soliton parameters 
/LtS;(0), !/fe(0), ^fc(O), i5fc(0) which characterize the region in the soliton parameter space 
responsible for the 7V-soliton bound states. 



3 The perturbed NLS and perturbed CTC 

We will consider several specific choices R'^^^u] of perturbations, p = 1,2,... in 
ijTjl. In the adiabatic approximation the dynamics of the soliton parameters can be 
determined by the system (see |TJ for N = 2 and OH] for N > 2): 



dX 
~dt 

dt " + "''^ ' dt 

where = /i^ + ivk and xjf^ — + D^k^. The right hand sides of Eqs 

l(Tn|l - ((TT|l are determined by i?^^^ [u] through: 



2^, + 4^', ^ = 2{^^l + .l) + xi^\ (11) 



1 r dzk 



M; 



r^/^^H^if !...(<[*-*). (13) 

^i"^4£^'''(<'i«i-'")' (") 

Inserting (^UJi lfTTl> h^to Q we derive: 

^gi = ^Av,Xk + -AA(^^ + 2*6 (M^i^ + ^N^^A 
at vq \ / 

+ z(2AoSi^)-4^)-A'(^)), (16) 



i=i i=i J=i 

In deriving eq. we have kept terms of the order /S.Vk — ©(y^) and neglected 
terms of the order O(eo)- The perturbations result in that and /ip niay become 
time-dependent. Indeed, from 1)10(1 we get: 
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The small parameter eo can be related to the initial distance ro = |^2 — fi|t=o 
between the two solitons. Assuming i>i^2 — we find: 



dx\uf{x,Q)uWxM ~ Si^oroe-^"'"''' . (18) 



In particular, p8|) means that eo — 0.01 for tq ~ 8 and z/q = 1/2. 

We assume that initially the solitons are ordered in such a way that ^k+i — Cfc — 
ro. One can check [HI that iV^^^ ~ Aljf'' ~ exp(-2i^o|fc - p\ro)- Therefore the 
interaction terms between the fc-th and fc± 1-st solitons will be of the order of e^^'^f'o ; 
the interactions between fc-th and fc ± 2-nd soliton will of the order of e-4i/oro ^ 

g-2i/oro 

The terms S^"'' , xj^^ are of the order of rg exp(— 2i'oro), where a = or 1. However 
they can be neglected as compared to Jlk and i^k, where 

Mfc = Mfe - A^o - V^o, I'fc = i^k - J^o ^ V^o, (19) 

The corrections to Njf"^ , . . . , coming from the terms linear in u depend only on 
the parameters of the fc-th soliton; i.e., they are 'local' in fc. The nonlinear terms in 
u present in produce also 'non-local' in fc terms in N^^\ .... 

3.1 Nonlinear gain and second order dispersion 

Consider the NLS eq. Q with 

R[u] ^ cou + C2Uxx + do\u\'^u, (20) 

where cq, C2 and do are real constants, see |12| . Another important factor is the order 
of magnitude of the perturbation coefficients Cq, C2 and do in H20|l . If we take them 
to be of the order of eo we find that the A^-soliton train evolves according to: 

^ = C/oo + le^^o^ (e'3'=+-Q'= - e^^-Q"-^) , (21) 

where for ^o = we get Uqo = — (^Scq + Si^qiIq — 4c2i'o). This form of perturbed 
CTC (I2U can be solved exactly: 

Qk{t) = ^Uoot'^ + Voot + Q'itHt), 

where Q'j^\t) is a solution of the unperturbed CTC and Vqq is an arbitrary constant. 
In this case the effect of the perturbation will be an overall motion of the center 
of mass of the A^-soliton train. The relative motion of the solitons will remain the 
same. For larger values of the coefficients co, C2 and do, e.g., of the order of y/eo the 
corresponding dynamical system is more complicated and has to be treated separately. 

3.2 Quadratic and periodic potentials 

Let iR[u] — V{x)u{x,t). Our first choice for V{x) is a quadratic one: 

V^^^x) ^V2X^ + V1X + V0. (22) 
Skipping the details we get the results: 

Ni^^=0, Af« ^-1/2^,-^, (23a) 
^i'^ = 0. d['^ = V2 f ^ - V^^k Vo, (23b) 
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and X^ip = D^jp . As a result the corresponding PCTC takes the form 

fetM = _4.o (eQ^.-Q^ - e^-e-.) _ v,^, - (24) 

^ = -4.o(M.+-.)-^i^i^^-^E<\ (25) 
If we now differentiate l|25(l and make use of H24(l we get [H]: 

It is reasonable to assume that V2 — O{eo/N); this ensures the possibility to have the 
A^-soliton train 'inside' the potential. It also means that both the exponential terms 
and the correction terms M^^'' are of the same order of magnitude. From eqs. (|24|l 
and (|25(l there follows that dvo/dt = and: 

where /io is the average velocity and — jj center of mass of the 

A^-soliton train. The system of equations for V2 > has a simple solution 

^ia{t) = ^ioo cos($(0), (o{t) = sin($(i)) - (27) 

where $(t) = \/2\^t + $0, and /ioo and $0 are constants of integration. Therefore 
the overall effect of such quadratic potential will be to induce a slow periodic motion 
of the train as a whole. 

Another important choice is the periodic potential 

V'-^\x) ^ Acos{nx + Qa), (28) 

where A, and flo are appropriately chosen constants. NLS equation with similar 
potentials appear in a natural way in the study of Bose-Einstein condensates, see [17) . 

For N = 2 solitons the corresponding Karpman-Solov'ev system was derived in 
PT] . For N > 2 we obtain the PCTC where the integrals for Nk, Mk, and Dk are 
equal to [S]: 

<' = o, „f. = ^_i^...,na, + n„, (2,) 

Hf^O, + (30, 

where Zk — 7r51/(4;^fc). These results allow one to derive the corresponding perturbed 
CTC models. Again we find that dv^/dt — 0. 



4 Analysis of the perturbed CTC and comparison 
with numerical simulations 

The dynamics of an individual soliton in a train is determined by the combined action 
of external potential and the influence of neighboring solitons. The interaction with 
neighboring solitons can be either repulsive, or attractive depending on the phase 
relations between them. Particularly, if their amplitudes are equal and the initial 
phase difference between neighboring solitons is tt (as considered below) they repel 
each other giving rise to expanding motion in the absence of an external field |21 IS| ■ 
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The external potential counterbalances the expansion, trying to confine solitons in 
the minima of the potential. It is the interplay of these two factors - the interaction of 
solitons and the action of the external potential, which gives rise to a rich dynamics 
of the iV-soliton train. 

To verify the adequacy of the perturbed CTC model for the description of the 
iV-soliton train dynamics in external potentials we performed comparison of predic- 
tions of corresponding perturbed CTC (PCTC) system and direct simulations of the 
underlying NLS equation Below we present results pertaining to a matter- wave 
soliton train in a confining (i) parabolic trap and in (ii) a periodic potential modelling 
an optical lattice. 

Here we present the numerical verification of the PCTC model. The perturbed 
NLS eq. ^ is solved by the operator splitting procedure using the fast Fourier 
transform |24| . In the course of time evolution we monitor the conservation of the 
norm and energy of the N-soliton train. The corresponding PCTC equations are 
solved by the Runge-Kutta scheme with the adaptive stepsize control [53] . 

The evolution of a 7V-soliton train in the absence of potential {V{x) = 0) is well 
known, see e.g. These papers propose a method to determine the asymptotic 

dynamical regime of the CTC for a given set of initial parameters J/fc(0), /Xfe(O), ^fe(O) 
and Sk{0). Below we will use mainly the following set of parameters: 



z/fe(0)-l/2, /ifc(0) = 0, 6+i(0)-6(0) = ro, 4(0) = fcTr, (31) 
MO) = 1/2, ^ik{0) = 0, a-+i(0)-a-(0) = ro, 4(0) = 0. (32) 



These two types of initial conditions (IC) are most widely used in numeric simulations. 

In the absence of potential the IC 1)3 Ifl ensure the so called free asymptotic regime, 
i.e. each soliton develops its own velocity and the distance between the neighboring 
solitons increases linearly in time. At the same time the center of mass of the soliton 
train stays at rest (see the left panel of Fig. Under the IC (|32|l the solitons 
attract each other going into collisions whenever the distance between them is not 
large enough. 

From mathematical point of view the IC (|31|l reduce the CTC into a standard (real) 
Toda chain for which the free asymptotic regime is the only possible asymptotical 
regime. On the contrary the IC H32|l lead to singular solutions for the CTC (see 
121 El 123)- The singularities of the exact solutions for the CTC coincide with the 
positions of the collisions. 

Below we will study the effects of the potentials for both types of IC. One may 
expect that the quadratic potential will prevent free asymptotic regime of IC (|31|l no 
matter how small V2 > is and would not be able to prevent collisions in the case 
of IC H32I) . The periodic potential, if strong enough should be able to stabilize and 
bring to bound states both types of IC. 

4.1 Quadratic potential 

For the quadratic external potential V{x) = V2X^ + V\x + Vb the perturbed CTC 
equations in terms of soliton parameters have the form: 




(33) 



(34) 



da 



(35) 



dt 
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dS ( ^ 



dt 



(36) 



where fik, i^k, and Sk for k ~ 1,. . . ,N are the AN sohton parameters, see eqs. 

©-Q- 

The effect of the quadratic potential on the A^-sohton train with parameters H31() is 
to balance the repulsive interaction between the solitons, so that they remain bounded 
by the potential, as illustrated in figures: the right panel of and in[21 The quadratic 
potentials are supposed to be weak, i.e. we choose V2 so that 

^2Cw(0) < ^0, V^2e?(0) < 1^0. (37) 
Figures n and m show good agreement between the PCTC model and the numerical 
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Figure 1: Left panel: 9-soliton train with initial parameters as in eq. (|31|l with tq = 8 
in the absence of a potential goes into free asymptotic regime. Solid lines: direct 
numerical simulation of the NLS equation Q|; dashed lines: £,k{t) as predicted by the 
CTC equations ((23)-® with Vq = Vi = V2 = and ro = 8. Right panel: Evolution 
of a 9-soliton train with the same initial parameters in the quadratic potential V{x) = 
V2X^ with Vi ~ 0.00005. Sohd hnes: direct numerical simulation of the NLS equation 
f^l; dashed lines: solution of the PCTC equations H33|) - (|36|) . Initially the train is 
placed symmetrically relative to the minimum of the potential at a; = 0. 



solution of the perturbed NLS equation (Q. They also show two types of effects of 
the quadratic potential on the motion of the iV-soliton train: (i) the train performs 
contracting and expanding oscillations if its center of mass coincides with the mini- 
mum of the potential, (ii) the train oscillates around the minimum of the potential 
as a whole if its center of mass is shifted. In the last case contracting and expanding 
motions of the soliton train is superimposed to the center of mass dynamics. As one 
can see from the figures the period of this motion matches very well the one predicted 
by formula H27() . Indeed, from eq. H27|l it follows that the period period of the center 
of mass motion is T = 27r/-v/2V2. For the parameters in fig. [2 we have T ~ 628 (for 
9-soliton train), T ~ 140 (for 3-soliton train). Similar is the dynamics also for the 
7-soliton train on Fig. 01 for the parameters choosen there we have T = 314.2, in good 
agreement with the numerical simulations. The direct simulations of the NLS equa- 
tion shows that stronger parabolic trap may cause merging of individual solitons 
at times of contraction, and restoring of the original configuration when the train 
is expanded. This behavior reminds the phenomenon of "missing solitons" observed 
in the experiment However, this situation is beyond the validity of the PCTC 
approach. 
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Figure 2: Harmonic oscillations of a N-soliton train initially shifted relative to the 
minimum of the quadratic potential V{x) = V2X^. Left panel: 9-soliton train, V2 — 
0.00005. Right panel: 3-soliton train, V2 — 0.001. The IC of the both trains is given 
by with ro = 8. In both panels solid lines correspond to direct simulations of 
the NLS equation QJ, and dashed lines to numerical solution of the PCTC equations 
- 



4.2 Periodic potential 

Another external potential in which the iV-soliton train exhibits interesting dynamics 
is the periodic potential of the form V{x) = A cos{ilx + i^o) . This case also may have 
a direct relevance to matter - wave soliton trains confined to optical lattices. The 
PCTC system in terms of soliton parameters has the form: 



dt 

dvk 
dt 



d^ 

dt 
d6k 

dt 



IGiyl (^e-^Mi.+i-i.) cos (2^^0(6+1 - 6) + 4 - 4+i) 
e-2.o(?.-«.-i) cos (2/io(a - 6-1) + Sk-i - 4)) + Mf 
161^0' (e-2''«(«''+i-«'=) sin (2Aio(&+i - 6) + Sk - 
g-2.o(a-C.-i) sin (2fioi^k - 6-1) + 4-1 - 4)) , 

2Mfc, 

2i^,l + 4) + D^^\,.k), 



(38) 
(39) 

(40) 
(41) 



where Ml >{i^k) and D)^'{uk) are given in ^ and 

Each soliton of the train experience confining force of the periodic potential and 
repulsive force of neighboring solitons. Therefore, equilibrium positions of solitons 
do not coincide with the minima of the periodic potential. Solitons placed initially 
at minima of the periodic potential (Fig. ^ perform small amplitude oscillations 
around these minima, provided that the strength of the potential is big enough to 
keep solitons confined. As opposed, the weak periodic potential is unable to confine 
solitons, and repulsive forces between neighboring solitons (at phase difference tt) 
induces unbounded expansion of the train. In the intermediate region, when the 
confining force of the periodic potential is comparable with the repulsive forces of 
neighboring solitons, interesting dynamics can be observed such as the expulsion 
of bordering solitons from the train, as shown in the left panel of Fig. El This 
phenomenon, revealing the complexity of the internal dynamics of the train, can be 
explained as follows. Each soliton performs nonlinear oscillations within individual 
potential wells under repulsive forces from neighboring solitons. When the amplitude 
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Figure 3: Dynamics of a 7-soliton train placed asymmetrically relative to the minimum 
of the trap Vix) = 0.0002a:^. Solid lines: results of direct numerical simulations of 
the NLS equation Dashed lines: result of solution of the PCTC system H33|l - 
l|36|l for the center of mass ^i. The parameters of solitons are the same as in H31|l with 
ro = 8. The initial shift of the soliton train relative to the minimum of the parabolic 
trap is IOtt. 
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Figure 4: Solitons (continuous line) remain confined around the minima of the periodic 
potential Vix) — j4cos(x) (dashed line) performing small amplitude oscillations if its 
strength is big enough A = 0.1. 



of oscillations of particular solitons grow and two solitons closely approach each other, 
a strong recoil momentum can cause the soliton to leave the train, overcoming barriers 
of the periodic potential. In Fig. Elthis happens with bordering solitons (the other 
solitons remain bounded under long time evolution). It is noteworthy to stress that 
this phenomenon is well described by the PCTC model, as is evident from Fig. |S1 left 
panel. 

On the right panel of the same figure we have similar IC as in (|31|l and we have 
choosen again the initial positions of the solitons to coincide with the minima of the 
periodic potential V{x) = ^cos(na: + flo); i-e. — 27r/ri. The values of ^ = 0.0005 
and ro = 9 in the right panel of Fig. [Slnow are such that the solitons form a bound 
state. Therefore for any given initial distance ro there is a critical value Acr(ro) for 
A such that for A > Acr{ro) the soliton train with IC (|31|l will form a bound state. 

In contrast to the quadratic potentials, the weak periodic potential is unable 
to confine solitons, and repulsive forces between neighboring solitons (at i'fe(O) = 
1/2, (5fe(0) = fcTr) induces unbounded expansion of the train similar to what was 
shown in the left panel of Fig. ^ 

The periodic potential can play stabilizing role also for the IC (|32|) . when the zero 
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Figure 5: Left panel: The expulsion of solitons from the train, as obtained from 
direct simulations of the NLS equation (solid lines), and as predicted by PCTC 
system H38|l - ()41|l for the center of mass (dashed lines). The IC of the 7-soliton 
train are given by lpTT|l with ro = 8; the parameters of the periodic potential V{x) — 
Acos(r2a; + fio) are A = 0.001, fl — ir/i, = 0. Right panel: Oscillations of the 
5-soliton train with IC given by ()31|l with tq = 9; in a moderately weak periodic 
potential, A — 0.0005, — 27r/9, JIq = 0. Solid and dashed lines correspond, 
respectively, to numerical solution of the NLS eq. ^ and PCTC system (|38f) - (|41ll . 



phase difference between neighboring solitons correspond to their mutual attraction. 
If the periodic potential is strong enough, solitons do not experience collision. The 
weak periodic potential cannot prevent solitons from collisions, which eventually leads 
to destruction of the soliton train, as illustrated in Fig. |S| Again for any given initial 
distance ro there will be a critical value A'^^{ro) for A such that for A > yl^^(ro) the 
soliton train with IC l|32|) will form a bound state avoiding collisions. 

4.3 Tilted periodic potential 

Now we consider the dynamics of a N - soliton train in a tilted periodic potential, 
which is the combination of periodic and linear potentials 

V{x) ^ Acos{nx + no) + Bx. (42) 

This potential is of particular interest in studies of Bose-Einstein condensates. A 
train of repulsive BEC loaded in such a potential (where the periodic potential was 
a ID optical lattice and the linear one was due to the gravitation) exhibited Bloch 
oscillations [22]. At each period of these oscillations condensate atoms residing in 
individual optical lattice cells coherently tunneled through the potential barriers. 
This was the first experimental demonstration of a pulsed atomic laser 27 . Recently 
a new model of a pulsed atomic laser was theoretically developed in _28^, where the 
solitons of attractive BEC were considered as carriers of coherent atomic pulses. 

Controlled manipulation with matter - wave solitons is important issue in these 
applications. Below we demonstrate that solitons of attractive BEC confined in tilted 
optical lattice can be flexibly manipulated by adjustment of the strength of the linear 
potential. In Fig. [3 we show the extraction of different number of solitons from the 
5-soliton train by increasing the strength of the linear potential B, as obtained from 
direct simulations of the NLS equation (1) and numerical integration of the PCTC 
system (38) - (41) with 
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Figure 6: Dynamics of a 5-soliton train with zero phase difference between neighboring 
sohtons, in the periodic potential V{x) = Acos(r2x + flo) with fl = 7r/4, flo = 0, and 
ro = 8. Left panel: When the periodic potential is strong enough A = 0.02, the 
N-soliton train remain confined, each soliton performing small amplitude oscillations 
around the minima of individual cells. Right panel: Weaker periodic potential A = 
0.01 cannot prevent solitons from collisions, which destroy the train. Solid and dashed 
lines correspond, respectively, to numerical solution of the NLS cq. and PCTC 
system §^ - (gH). 
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As is evident from Fig. |7| the PCTC model provides adequate description of the 
dynamics of a N-soliton train in a tilted periodic potential. A small divergence be- 
tween predictions for the trajectory of the left border soliton in Fig. \7\ (d), is due 
to the imperfect absorption of solitons from the right end of the integration domain. 
Reflected waves enter the integration domain and interact with solitons, which causes 
the discrepancy. 



5 Conclusions 

We have studied the dynamics of the A-soliton train confined to external fields 
(quadratic and periodic potentials). Both the analytical treatment in the framework 
of the PCTC model, and numerical analysis by direct simulations of the underlying 
NLS equation show that the PCTC is adequate for description of the A-soliton inter- 
actions in external potentials. Relevance of this study to the research on matter- wave 
soliton trains in magnetic traps and optical lattices is briefly mentioned. 

Like any other model, the predictions of the CTC should be compared with the 
numerical solutions of the corresponding NLEE. Such comparison between the CTC 
and the NLS has been done thoroughly in ^,'3','5^ and excellent match has been found 
for all dynamical regimes. This means that the CTC may be viewed as an universal 
model for the adiabatic A-soliton interactions for several types of NLS. 

For the perturbed CTC equations such comparison has been just started; the 
good agreement shown in the figures above supports the hope that the region of 
applicability of PCTC can be widened. 

More detailed investigation of the A-soliton train interactions under different types 
of external potentials and for diffetent types of initial soliton parameters will be 
published in subsequent papers. 
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Figure 7: Controlled withdrawal of solitons from the 5-soliton train by adjusting 
the strength of the linear potential 142|l with parameters: A = 0.0005, Vl — 27r/9, 
Oo = 0. Depending on the tilt, different number of solitons can be pulled out of 
the train: (a) one soliton at B = -0.00003, (b) two sA. B = -0.00011, (c) four at 
B = —0.0002 and (d) five at B = —0.0003. The initial phase difference and separation 
between neighboring solitons in the train are, respectively, tt and 9. Initially the 
train is shifted by — IOtt with respect to a: = for graphical convenience. Solid and 
dashed lines correspond, respectively, to direct simulations of the NLS equation 
and numerical integration of the PCTC system (|38|l - H41|) . 
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